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OF THE SCHWINGER-DYSON EQUATIONS 
FOR THE FERMION AND PHOTON PROPAGATORS 



' A. Kizilersii^'^ and M.R. Pennington^ 

O 

^ I ^ Special Researcli Centre for tlie Subatomic Structure of Matter, 



p ^. Department of Physics, 

Adelaide University 
Adelaide, 5005, Australia 



< 
(N 



> 
cn 
oo 

m 
o 



2 



Institute for Particle Physics Phenomenology, 



JlLi' Physics Department 

^ . Durham University 



Durham DHl 3LE, U.K. 
ABSTRACT 

In principle, calculation of a full Green's function in any field theory requires knowledge of 
the infinite set of multi-point Green's functions, unless one can find some way of truncating 
^ . the corresponding Schwinger-Dyson equations. For the fermion and boson propagators in 

QED this requires an ansatz for the full three point vertex. Here we illustrate how the 
' properties of gauge invariance, gauge covariance and multiplicative renormalizability impose 
;h ' severe constraints on this fermion-boson interaction, allowing a consistent truncation of 



the propagator equations. We demonstrate how these conditions imply that the 3-point 
vertex in the propagator equations is largely determined by the behaviour of the fermion 
propagator itself and not by knowledge of the many higher point functions. We give an 
explicit form for the fermion-photon vertex, which in the fermion and photon propagator 
fulfills these constraints to all orders in leading logarithms for massless QED, and accords 
with the weak coupling limit in perturbation theory at 0{a). This provides the first attempt 
to deduce non-perturbative Feynman rules for strong physics calculations of propagators in 
massless QED that ensures a more consistent truncation of the 2-point Schwinger-Dyson 
equations. The generalisation to next-to-leading order and masses will be described in a 
longer publication. 
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I. INTRODUCTION 



Solution of the Schwinger-Dyson equations (SDE) for any field theory would constitute the 
complete determination of that theory and every possible measurable quantity would be 
known. Even though it is nearly 60 years since these field equations were first derived 
S, Q, 0], we are far from obtaining their solution even for a relatively simple theory like QED. 
Progress has been hampered by the very structure that makes field theory interesting, namely 
that the Schwinger-Dyson equations form an infinite nested set. Each n-point function 
must be multiplicatively renormalizable and, in a gauge theory, respects gauge invariance. 
To achieve this, the solution even for the 2-point functions (the propagators) appears to 
require knowledge of all the other n-point functions. Consequently, studies in gauge theories 
have resorted foremostly to a perturbative approximation, in which each Green's function 
is expanded to a given order in the coupling squared. Or as an approximation to non- 
perturbative physics, simple (even simplistic) ansatz have been used for the 3-point function 
to allow the fermion propagator to be investigated. In return dynamical mass generation 
has been studied in the rainbow approximation 

aaii 

[lol . [ill . [12I and some level of 

understanding of when chiral symmetry breaking can occur has been reached. While valuable 
for gaining intuition, this is no substitute for a genuine non-perturbative study. While formal 
results on gauge invariance and multiplicative renormalizability (MR) have long been known 
using the gauge technique of Salam, Delbourgo , [isl and others, this method has 

not proved useful for providing equations that can be readily solved either analytically 
or numerically. Here, an alternative approach, an attempt to develop non-perturbative 
Feynman rules, has proved more fruitful. The aim is to write down explicit representations 
for the effective n-point functions, in particular, for the 3-point function, which ensures 
that the solutions of the Schwinger-Dyson equations for the 2-point functions respect gauge 
invariance and are multiplicatively renormalizable 16l. Il7|. 



What has previously impeded the practical study of the Schwinger-Dyson equations has been 
the need to handle overlapping divergences that dramatically complicate the renormalization 
of the equations. The present approach overcomes this difficulty by requiring that the 2- 
point functions must be multiplicatively renormalizable and no overlapping divergences can 
thereby occur. This procedure is genuinely non-perturbative and is not readily relatab 
att emp ts at summing subsets of Feynman graphs with these same properties 
H. 
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The first of such non-perturbative studies has been in the case of quenched QED 
^1— 11— 11— 11— 11— I 1—11—11—11—11— 



25 



26,^ 



28 



29 



32 



33 



34 



12 



23 



35| — that is QED in which the explicit factor of 



Np multiplying the fermion loop corrections to the photon propagator is set equal to zero. 
Then a form for the fermion-boson vertex that satisfies the Ward identity, the Ward-Green- 
Takahashi identity 36|] and renders the fermion propagator multiplicatively renormalizable, 
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has been written down explicitly [18|, |37j. While the form is non-perturbative, the fact 
that it must agree with perturbation theory in the weak coupling regime is a key pointer 
to the ultraviolet structure, expressed in terms of logarithms of momenta. The purpose 
of the present paper is to extend this study by developing the constraints that have to be 
fulfilled in the case of massless unquenched QED to ensure both the fermion and photon 
propagators are multiplicatively renormalizable (at least as far as leading logarithms are 
concerned) . 

In general, the full fermion-boson vertex has 12 components, all of which are in principle 
independent, though one is forced to be zero by gauge invariance. The fermion and photon 
propagators do not require complete knowledge of the full complexity of this structure, but 
just 2 projections that arise in the Schwinger- Dyson equations for these 2-point functions. 
We present a simple solution to the constraints from multiplicative renormalizability. While 
the general structure of the full vertex is not complete, the projections within the SDEs for 
the 2-point functions have no freedom. 

While it is clear that the full 3-point function must involve knowledge of the 4-point kernel 
and higher-point functions, as far as its role in the equations for the propagators is concerned. 



this is not the case. Thus it can be that the effective 3-point function involves on 
point functions. A clue to this is provided by the Ward-Green- Takahashi 



38 



39 



y the full 2- 



40| identity, 



which tells us that part of the 3-point vertex (often called the longitudinal part) is precisely 
fixed by the fermion propagator alone. Moreover, a hint that the remaining transverse 
part may be similarly constrained is the fact that the vertex and fermion wavefunction 
renormalization have common renormalization factors (Zi = Z2) as a consequence of gauge 
invariance. Thus the transverse part must know about the fermion propagator functions 
too. How, this works in full QED is what we investigate in this paper. 

In Sect. 2 we consider the structure of the fermion-boson vertex and its ultraviolet be- 
haviour. In Sect. 3 we compute the Schwinger- Dyson equations for the fermion and boson 
propagators. In Sect. 4 we deduce the ultraviolet structure imposed by multiplicative renor- 
malizability. Sect. 5 gives the constraints on the vertex imposed by MR conditions. The 
pattern of constraints indicates a general analytic form for the transverse part of the vertex 
structure. In Sect. 6 we deduce a solution to these constraints involving the full fermion 
wavefunction renormalization. The vertex in the weak coupling limit is studied in Sect. 7 
and the restrictions it imposes derived. In Sect. 8 we conclude and outline a programme for 
future work. Since this procedure is rather complicated, we show in Fig. 1 a flow diagram 
of this calculation. 



3 



Non-perturbative massless QED 

t 



Schwinger-Dyson equations 



Truncation is needed 
Make an ansatz for 3-point vertex 



T 



As a result of Gauge 
Invariance Ward Identities 
must be fulfilled 



basis tensors are 
given by Ball-Chiu 



Using charge conjugation, 

other information and. 
perturbative expansion of 
Ti are suggested 



The ansatz must satisfy criteria, 
which the full vertex itself 
satisfies 




Full vertex is divided into 
longitudinal and transverse 

parts 





Longitudinal verteBTls fixed 

by Ball-Chiu [36| 



^ BC 



Transverse part left to be 
determined. The vector 
structure of the vertex leads to 



i J, 



Ti, the coefficient functions 
are the only unknowns 



Substitute this vertex into 
the coupled SD-equation 



Calculate the general form of 

multiplicatively 
renormalisable F and G 



Perturbative calculation and 
all other vertex information 



Solve SD-eqn. for the 1/F and 1/G 
in terms of the constants of 
the perturbative expansion of Ti 



Find the constraints on the vertex 

function imposed by 
multiplicative renormalisation 

and generalize these constraints 



Construct non-perturbative vertex ansatz 



FIG. 1: Flow diagram of the Schwinger-Dyson calculation presented here. 
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II. VERTEX AND PROPAGATORS AND THEIR RENORMALIZATIONS 
DEFINED 



The two key constraints on the fermion-boson vertex are provided by the gauge invariance 
of the theory and by muhiphcative renormahzabihty. Here we begin with the first of these 



and describe the importance of the Ward- Green- Takahashi identity |38|, |39|, |40|. Though 



this is well known, it forms the essential background allowing us to establish our notation. 




FIG. 2: Fermion-boson vertex. 

The vertex, displayed in Fig. [21 is a function of the two independent momenta flowing 
through the vertex. We take these to be the fermion momenta, k and p. The vertex function 
is r'^{k,p; q) with q = k — p . It is well-known that the coupling of two spin-1/2 particles 
to a spin-1 boson involves 12 independent vectors, of these eight are transverse to the boson 
momentum q. The structure of the four (longitudinal) components are constrained by the 
Ward-Green-Takahashi identity (WGTI) 

q^T,{k,p; q) = S^\k)-S^\p) , (1) 

where Sf{p) is the full fermion propagator carrying momentum p. In general 

1 

where F{p'^) (or Aip^) = 1/F{p'^)) is the fermion wavefunction renormalization and M.{p'^) 
(or BijP') = M ijP') / F ijp')) is its mass function. The bare fermion propagator is just 
S^{p) = 1/ (^ — mo). From the form of this propagator, the Ward-Green-Takahashi identity, 
Eq. ([T]), contains terms with both one and no gamma matrices, so that the vertex component 
involving two through a'^'^ = |[7^,7''] must be zero. Thus in a gauge theory there are in 
fact 11 independent non-zero vectors in terms of which to decompose r^(p, k]q). Of these. 
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six occur if the fermions are massless as we consider here, i.e. A4{p'^) = 0. Eq. ([T|) has a 
well known zero photon momentum limit; the Ward identity : 

T,{p,p; 0) = hm r,(p, k;q) = ^^^f^ . (3) 
fc >P op^ 

The full vertex can be divided into longitudinal and transverse components 

T>^{p,k;q) = Tiip,k;q)+Tf;,ip,k;q) , (4) 

where 

q^Tt^{p,k;q) = . (5) 

We demand that the longitudinal part alone is responsible for the vertex satisfying both 
Eqs. f iniSl) . This means that each component must be separately free of kinematic singular- 
ities, so that 

r^(p,p;0)=0 . (6) 



The longitudinal part is then defined, following Ball-Chiu 36|], to be 

^1{p,k,q) = r^^(p,A:,g), 



J2HP^k',q')mp,k;q) , (7) 



where 



^i(^'^^^^')=^:^ + :^J , Lnp,k;q)=^^, 
^^^^''^''^')=^(P^(f(W"^) ' Li^{p,k;q) = {k'^+pn{H 



i=l 



{k^ - p2) \^ F(fc2) 

A4(p^ q')=0, L>i{p, k- q) = a^'' {K + Pu) • (8) 

Crucially because of gauge invariance, this longitudinal component of the vertex is wholly 
determined by the fermion propagator. Moreover, it is this longitudinal component that 



gives the dominant ultraviolet behaviour of the vertex 41 1. 



Quite generally, the transverse vertex can be decomposed in the massless fermion case in 
terms of the remaining four basis vectors as : 

r^(p,A;;g)= n{p\k',q')Tt{p,k;q) , (9) 

1=2,3,6,8 
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where the are coefficient functions depending on momenta fc^,j9^ and g^, which are as 
yet undetermined, and the Tj are the basis tensors defined by Ball and Chiu 36|] — the 
modification of this basis by Kizilersii et al. 4l| does not affect these four vectors : 

T^{P. k; q) = {p^{k ■ q) - k^{p ■ g) ) + ^) , 

n{p,k-q) =q'Y-q'^, 

T^{p,k-q) =^'^(p^-e) + {p + kr^, 

Ts^{p,k;q) =-rkVaxu + ky-p^l^. (10) 

With these basis vectors, the Ti{i = 2, 3, 6, 8) are individually free of kinematic singularities 
at 0{a) in perturbation theory in any covariant gauge as shown in Ref. 4]|. It is these r^'s 
that are constrained by multiplicative renormalizability [3]. It is our key presumption that 
this will force these transverse components (or at least their projections in the Schwinger- 
Dyson equations for the 2-point functions) to depend only on propagator functions just like 
the longitudinal part of Eqs. ([71 [8]). 

What we can say about these coefficients? 

Here we discuss the fundamental constraints on the transverse vertex that follow from (z) di- 
mensional analysis, (ii) symmetry properties, {in) order of perturbation theory, (iv) gauge 
dependence and (v) renormalization : 

(z) The transverse vertex is dimensionless. Knowing the dimensions of the basis vectors from 
Eq. ffTOj) tells us the dimensions of the r^'s. With d = momentum^, then: 

1 



dim. of T2 



dim. of TJ^ 



24 



dim. of Tg 
dim. of 

[ii) The C-parity operation 
r2{e,p',q') = 

T,{k\p\q')= r,{p\k\q'\ 
r,{k\p\q')= -T,{p\k\q' 
T,{k\p\q')= n{p\k\q'\ 



d' 
d 
d 
d 



dim. of T2 
dim. of T3 
dim. of Tq 
dim. of Tg 



d^ 
1 

d ' 
1 

d ' 
1 

d ■ 



(11) 



42| on Eqs. (j?! ^ requires 

T,{p',e,q') , Xi{e,p',q') 

, X2{k\p',q') 

, Xs{e,p',q') 

, Hk\p\q^) 



\i{p\k\q') 
\2{p\k\q') 
Hp\k\q') 
-Hp\k\q^) 



(12) 
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[iii) At zeroth order in perturbation theory the full vertex is 7'^. Since at this order F = 1, 
we see from Eqs. ([71 [H]) that = 7^^, consequently, = 0. Thus the = 0{a) in 
perturbation theory. 

(iv) The propagator for the photon carrying momentum q is 



fJ.U 



+ e 



q' 



(13) 



where G{q^) is the photon renormalization function, ^ is the covariant gauge parameter and 
the A^^, is the transverse part of the photon propagator. The bare photon propagator, A 



1 



has = 1 in Eq. ^ 



Gauge covariance is expressed through the Landau-Khalatnikov-Fradkin (LKF) transforma- 



tions 



43 



44( 1 . These mean that once a Green's function is known in some gauge, then its form 



in all other gauges is determined. In general, this is, of course, only useful if we know the 
relevant Green's function precisely in some gauge. Nevertheless, the LKF transformations 
provides two key results we shall use. The first concerns the fermion wavefunction renormal- 
ization, F{p^), which can only depend on the covariant gauge through a unique factor of ^ in 
its anomalous dimension. The second fact is that the photon wavefunction renormalization, 
G{q^), must be gauge independent. Both of these requirements place restrictions on the 
form of the non-perturbative interactions. 



{v) In QED the ful 
well known 



propagators and the vertex function are all divergent. However, as is 
45| , one can define finite (renormalized) propagators and vertex function 
by absorbing these divergences into functions, Zi {i = 1,2,3). As usual we introduce field 
renormalizations : 



(14) 



where the subscripts R and denote renormalized and bare quantities, respectively. The 
latter are conveniently made finite by introducing an ultraviolet momentum cutoff A and 
the former renormalized quantities depend on the momentum scale /.i at which we choose 
to renormalize. The divergence of the fermion propagator is absorbed into Z2, the fermion 
renormalization function, by : 



(15) 



and similarly for the photon : 

A^,(p,/i) = Zi\f,,A)A%ip,A). (16) 

The gauge covariance of the photon propagator requires that the covariant gauge parameter 
is similarly renormalized : 

Cn = Z,-'C- (17) 

The divergence of the vertex function is cancelled by the factor Zi : 

r>,/x)= Zi(/i,A)rO(p,A), (18) 

with the above definitions, the coupling constant is renormalized according to, 

Z2 



7 v^se. (19) 
^1 



Making use of the Ward-Green- Takahashi identity [38|, l39|, |40] : 

Zi = Z2, (20) 

the coupling constant renormalization becomes 

e^ = Z.^/\. (21) 

As usual, we define a = e^/(47r), where ao, an denote the bare and renormalized couplings 
related to e and e^, respectively. 

What we want to determine are the constraints these renormalizations of the fermion and 
photon propagators impose on the transverse part of the fermion-boson vertex. The renor- 
malization of the 3-point vertex is proportional to fermion renormalization constant Z2^- 



This can be seen already in the longitudinal vertex from the WGTI (36|]. Consequently, the 
non-perturbative structure of the transverse component, and hence the Tj's, must be pro- 
portional to the inverse of the fermion wavefunction renormalization, i.e. Ti {F,G) ~ 
just as the longitudinal Aj's of Eq. (8) are. 

To go further, the basic idea is easily explained by considering the fermion propagator 
in quenched massless QED. The non-perturbative quantity is the fermion wavefunction 
renormalization F(p^, A^). Let us imagine expanding this perturbatively and just keeping 
leading logarithms, so that we have 

2 2 2 

Fip',A') = 1 + aoA,\n^ + alA2\n'^ + alAsln'^ + ■■■ , (22) 



then inserting such a form in the loop integral of Fig. [31 For this to be a solution of the 
Schwinger-Dyson equation, the equation has to deliver -F(p^, A^) with the same perturbative 
expansion as output. However, to be multiplicatively renormalizable, the coefficients An 
cannot be independent, but related by A2 = As = Af/6 and finally An = A^ /n\. This 

requirement places a severe constraint on the fermion-boson vertex. Since its longitudinal 
part is known, it is its transverse components that are constrained. The aim of this paper is 
to determine these conditions on the Xj of Eq. ([9]) for full QED. In general, these rj(p^, /c^, q^) 
functions can be written as a sum of terms, each with the correct dimensions, Eq. (fTT|) . 
symmetry properties, Eq. f[T21) . and renormalization requirements, as : 

n{p\ k\ q') = J2 h {P\ k\ q') T? (F, G) . (23) 

j 

Each of these r/s has been divided into two parts : a kinematic part encoded in fij, giving 
the right dimensions, Eq. (fTTj) . which depends on momenta squared, and a functional part, 
rP, that is assumed only to know about the fermion and photon renormalization functions 
F and G at k^, p^ or q^. Such a form would provide a genuine non-perturbative construction, 

rf k\ q') = J2 [fT' k\ q') rr''^'^ (F, G) + f^^ {p\ k\ q') r^^^^ (F, G) , 
j 

Tr'\p\ k\ q') = (P'' k', q') rr'^'^ [F, G) + f^^ {p\ k\ q') rr^^^ [F, G) . 

(24) 

The forms of the Xj's are structured such that the integrals are soluble. First we deal with 
the kinematic factors for each r^'s, which are included in the following way: 
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+ 



(A;4 - p4) 
2 



1 

1 

1 

" (A;2+p2) 
(A;2-p2) 



/52 +72 
^2 + £2 



2k-p 

J^2 _|_ p2 



/^S + 73 

^3 + es 



A;2 + p2 
2k-p 



^2 _|_ p2 



{k^+p^y 
1 

(/i;2 — 
1 



/^e + 76 
^6 + ee 



^2 _|_ p2 

2k-p 



A;2 + p2 
2A;-p 



/^S + 78 
^8 + 68 



/c2 + p2 

2k-p 



^2 _|_ p2 

2A;-p 



A;2 + p2 



Tr'\p\k\q^) 
Tr{p\k\q') 

Tt\p\k\q^) 
rriP^k^q') 

Tt\p\k\q') 
Tr{p\k\q') 



(25) 



The factor 2 in the numerator of T2 is merely for later convenience and superscript "M" 
stands for Minkowski space. The kinematic factors, f^f"^' play two roles: first to ensure 
that each of t^^"^' is dimensionless, and to define the appropriate symmetry of these 
functions under the interchange of k^p. To make the problem tractable we do not include 
g2 dependence in the denominator factors. However, the dimensions and symmetry of the 
^sym, anti .^^ course, maintained by multiplying by a factor of q^ jik^ + p^\ Such a factor 
can be rewritten as \ — 2k-pj {k'^+p'^), and this is the origin of the inclusion of the (3i, 7^, 6i, 
terms in Eq. (1251) . 

The rf"** and r*^™ are antisymmetric and symmetric under k"^ ^ p"^, respectively. The 
^sym, anti assumed to be solely functions of the fermion and boson renormalization func- 
tions F and G, with consequently simplified dependence on k'^, p^ and q^. Since here we 
expand these functions in terms of leading logarithms, it is helpful to note that combinations 
like log(A;2/p2) antisymmetric, while \og{q'^ /k'^p^) is clearly symmetric under the inter- 
change of k and p, with each power of a "log" being multiplied by a factor of uq. Such forms 
are the basis for the leading logarithmic expansion of the r/^*"' Before renormalization, 
these will depend on the ultraviolet cutoff A, and we can represent the r*^™' by: 
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oo oo 



m=l n,r=0 



A2 



ao In 



A2 



A2 



fc2 p2 ^ 



(26) 



oo oo 



m=0 n,r=0 



ao In 



p 
A2 



2\ m- 



-,2 \ ^ 



tto In 



A2 



A;2 



2 \ »■ 



(27) 



The fact mentioned earlier that the zeroth order vertex contribution comes from the lon- 
gitudinal component, j'^, imposes the condition that there can be no leading order term 
in any transverse component. Consequently Sqqqq = 0. It is important to note that the 
coefficients A and S are constants in the above expressions and these depend on indices 
m, n, r. These are labelled by mnrr to make it easy to read off that such a term contributes 
at C(a^+"+"+"). Expanding Eqs. ([23 EZD to C(a|]): 

[p ,k ,q ) =ao Ai 



uooo 



k J) 

A2 A2 



On 



+ al< 



A 



2000 



fc2 



in^ - In^ i;^ 



2^ 



A2 



A2 

2 



A;2 p2 
+^iioo(ln^-ln-^ jln^ 



^3000 ( ^ - In^ ^ 



+^2100 In 



"0 



1200 



>1000 



fc2 



1^X2 "^''a2)^'''a2+-^1°^1 



A2 

.2 



A^ r" A2 



In —J - In — J 
A2 A2 



A;2 p2 



, A:2 p2 



+2 5oioo In 



(28) 



>200ol In ^+ln 



+2 5o2oo 

In 



2SL 
A2 



+ ao< 



-5iioo ( ^2 X2 ' 



A2 



53000 (^In^^+ln^l^j 

/ j^2 p2 \ ^2 

+52100 ( In'^+ln^— j In ^ 



+5ioii^ln-^ +ln|2 



A:2 p2 
+2 5ooiiln^In^ 

A2 



+2 5o3oo In'^ 



^2 

A2 



p. 
A2 



+ 5i2oo ( In + In ^ ) In^ -p 



P 



j^2 ^ p2 ^2 



In In -pj+2 5oiii In -pj In—^ In ■ 
A2 A2 A2 A2 A2 



+ 0{at) . 
(29) 
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One should keep in mind in the rest of this section that the sum of m, n, r, r adds up to the 
order of the expansion. Thus, for example at C(«o) one only has coefficients (^2000, ^1100) 
in rf"** and (.Saooo, '^0200, '^noo, ^oon) in rf^"". In turn, the dependence of A 
on ^ and Np can only happen such that the maximum power of each of them is m + n + 2r, 
i.e. the order of aQ too. 

As mentioned earlier the dominant ultraviolet behaviour of the vertex to O{ao) is given by 



the longitudinal component [4l|, Eq. ([71), the transverse vertex has no leading logarithms, 
i.e. (oq In" A^) terms must vanish. Consequently, in Eqs. (1271 [29]) the relation : 



at O{ao InA^ 



'^1000 + '^0100 — • (30) 



at 0{al In' A 



2 \2\ . 



'^2000 + '^0200 + '^0011 + '^1100 ~ • (31) 



at 0{al In^A^) : 



'^2100 + '^3000 + '^0300 + '^1011 + '^0111 + '^1200 ~ ^ • (^2) 



and in general at 0{aQ ln"A^) : 

u 

'^rn=[u~n~-2r]nrr — ) (33) 

nr=0 

must hold. 

Our aim is to determine the conditions on the constants Al^^rr iS^„rr 2 = 2, 3, 6, 8 im- 
posed by the fact that the fermion and photon propagators satisfy the appropriate Schwinger- 
Dyson equations and that these must be multiplicatively renormalizable. These constraints 
must, of course, be fulfilled by the full 3-point vertex. In the weak coupling limit, perturba- 
tive calculation of the relevant Feynman graphs will give explicit values for these constants. 
However, the Xj's that enter here determine not the full vertex, but projections defined by 
the Schwinger-Dyson equations of the next section. 
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III. UNQUENCHED SCHWINGER-DYSON CALCULATIONS 
A. Fermion propagator 




FIG. 3: Unquenched Schwinger-Dyson equation for fermion propagator. 

The Schwinger-Dyson equation for the fermion propagator displayed in Fig. [3] can be written 
as : 

- ^SF\p) = -^S'F\P) - [ 7^(-zer'^(p,fc;g)) ^SF{k) {-leY) ■ (34) 



Substituting the form of the longitudinal part of the photon propagator from Eq. (fT3|) and 
using the Ward-Green- Takahashi identity of Eq. ([1]), we can rewrite Eq. (134|1 as : 

iS-Ap) = ^^r(p)-e^^^^| T^{p.k-q)Sp{k)Y^lM 



^Sf\p) - £ ^1 r^(p, k; q) S,{k) Y 

+ SAp)SAk)l^y (35) 



The second term in the integrand being an odd integral gives zero : 



(2 7r)4 q 



4 



0, (36) 



if a translation invariant regularization is employed 2J]. After substituting the fermion and 



photon propagators, Eqs. ([21 [13]), explicitly in Eq. f[35l) . we obtain 
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(2vr) 



7 9 I fi'/l!^ 



(37) 



Multiplying this equation byp74, taking its trace and rearranging, we arrive at the following 
equation for the fermion wavefunction renormalization : 



F(p2,A2) 



1 + 



Tr ^{ T^{p,k-q) 1^^ F{k^) G{q^) [ g,. 



(38) 



We see this equation involves a particular projection of the full vertex T^. To make this 
explicit we substitute into this equation the general form given by the Ball-Chiu longitudinal 
part, Eq. ([7]), and the transverse component, Eq. ^ : 



1 



F{p^,A^) 



1 + 



ze 



d^k 

V(27r)" Jm W 



F{k' 



x< -^^Tr(^^^#) 

+ X^{p^k^q')G{q')^T 
+ A2(p^fc^g2)G(g2)Tr 
+ G(g')Tr 



fiL'^{p,k,q) fiY 9t,u 



^L!^{p,k,q) g^^ 



^Tt^{p,k,q) f^Y 



q2 



,(39) 



where d^k = 2n k^ dk"^ siv? ^ d"^ and ^ is the angle between the 4- vectors k and p. To 
perform these integrals, we move to Euclidean space using the Wick rotation [k^ — > ik^, 
ki — ki). After performing an explicit trace algebra in Eq. fl39|) and inserting the transverse 
vertex, F^, Eqs. ( iQllTOl) . with its undetermined r^'s, we obtain : 
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X 




-,{k'-k-p) 



+ F 




+ r^{p\k\q'){2^'^M'k-p] 
+ Ti{p\k\q'){-?.{k'-p')k.p] 

+ rs>^fc^g2){-2A2} ]\ , 



(40) 



where = [k ■ pY — k'^p^. 

Since multiplicative renormalizabilty is closely related to the ultraviolet behaviour of the 
Green's functions, we make a general perturbative expansion of the non-perturbative fermion 
and photon wavefunction renormalizations in powers of leading logarithms as follows : 



In this paper we will consider leading logarithms only in order to present the ideas and 
techniques and postpone to a future paper the more involved next-to-leading order. Of 
course, in perturbation theory the coefficients Bu have definite values. However, it is 
the general structure that multiplicative renormalizability determines. We substitute these 
expansions into Eq. (140|) in order to calculate this. The photon wavefunction renormalization 
G{q^) depend on the momentum q^ = k^ + p"^ — 2k ■ p therefore it has both a radial 
and an angular component. However, the angular dependent part of this quantity only 
contributes to 1/F{p'^) beyond the leading order, and so here we can simply approximate 
G{q^) with Giji^). We can then carry out the angular integration in Eq. (HOl) after inserting 




(41) 



M = 




OO 



(42) 



u=0 
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the coefficients of the basis tensors, i.e. Aj's and r^'s from Eqs. (18| [25|l : 



"0^ 



F(p2, A2) 



4 vr ./p2 



fc2 F(p2 



where 



3 ao 



2 Vi^(^^) 



(43) 



— anti 



p ^anti \ IP ^, \ ^anti i fp i ^, \ ^anti p 
P2 ^2 + (P3 - 73) ^3 + (P6 + 76) ^6 - P8 



. anti 



T 



sym 
f 



52 r^'"^ + (^3 - ^3) r^y^- + (^6 + e,) - 6s r, 



sym 



(44) 



To evaluate this expression, we have to insert the coefficients of the basis tensors, i.e. the 



anti, sym 



from Eqs. fl26l 1271) into Eq. fl43l) . A is the uhraviolet cut-off for the momentum k 



introduced in Eq. (l43l) in accord with Eqs. (fT5t[T6l[T8| 14111421) . One observes from Eq. (143|) 
that there is no contribution to l/F{p'^.,A^) from the Ai part of the longitudinal vertex, 
Eq. (I7j), but only from A2. On laboriously integrating Eq. (H3l) and using Eqs. (HTl H2l) we 
arrive at : 



F(p2,A2 



1 



(«+l) 



A2 



16 TT 



00 00 



u=0 t=0 



(u + t + 1) A2 



t=0 



(« + l) 



00 00 00 



Stt 



?i+t+m+n+2r+l |^«+t+m+n+2r+l P 



u=0 t=0 



m=l n=0 r=0 



A2 



"^mnrr 



(u + t + m + n + r + l) (n + t + ra + r + 1) 



00 00 00 



M+t+m+n+2r+l j^?i+t+m+n+2r+l P 



u=0 t=0 



n=0 m=0 r=0 



A2 



X s: 



+ 



(u + t + m + n + r + l) (u + t + n + r + l) 



where 



(45) 



•^mnrr — /^2 •^mnrr ~^ (/^3 73) •^mnrr ~^ ~^ 76) -^mnrr -^ri 



mnrr 



62 ^mnrr + ('^3 — ^s) '^^nrr + ('^6 + ^mnrr ~ ^8 ^mnrr 



(46) 
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In order to rearrange the infinite sums in Eq. (H5|) in terms of powers of Oq, we convert some 
of the infinite sums to finite sums : 



e 1 
4 7r F(p2^ 



(u + 1 



M=0 



A2 



3 

167 



oo 2 f " 1 

«=! ^^ a=l ^ ' 



u u—b 



b=l a=0 



u=l 



(47) 



where 

u b b—c c ( ^ ^ \ 



"y^-u.,— — ■ 

(48) 



u b b — c c ( ^ 

b=0 c=0 d=Oa=0 \^y2\^^")^^\ 



[^(-+^)-«+i]r'^-^^^' (49) 



with 



1 if j is even 
if j is odd . 

(50) 



The above expression for the fermion wavefunction renormalization, 1/F{jp,h?'), is the 
exact non-pertubative calculation for the massless fermions in a general covariant gauge at 
leading logarithmic order. In this equation the ^^„^^'s and iS^„^^'s are the constants to 
be constrained by multiplicative renormalization. For the purpose of explaining how this 
works, we will first implement it order-by-order then we generalize. To do this, we expand 
the fermion wavefunction renormalization, Eq. fHTj) . in O(tt^) : 
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1 + - (Xo^ In^ 

4:71 ] 



2, 2? 

an In 



A2 



e 3 
- + - 
2 8 



3 f-rf 

4 \ 1000 



5 



1000 



3 



1 A, 



+ -^{Ai + Bi) {a[qqq - 5/000 



A. 



2000 



^ -^llOO 



~7 i-' 



2000 



A '-'0200 



0011 



4i 4P 



19 1 

- A? 5i + —A2B1 Al B2 

4 ^ 16 16 



+ ^(A2 + i?2) U/000 - ^5 



/ 

1000 



-AiBi (a(qqq 



S 



1000 



{Al + 5i) (+ -^/ooo + g A^ioo 



1-^/ _ 1 0/ 
^'-'2000 g '-'0011 



0200 



1-r/ , 3-j-/ 9-j-/ 



+ ^ ■^lOll + g -^1200 + g -^2100 



A. 



3000 



, 1 

+ o '^0111 



3^/ _ 1^/ 

o '-'0300 



'-'1011 ^ ^'-'1200 ^'-'3000 



(51) 



Eqs. fl30H32l) have been input to obtain this expression. Eqs. fHSj Wfi [511) illustrate how the 
fermion wavefunction renormalization depends on the explicit form of the full 3-point vertex. 
As we shall see in Sect. IV, the expansion to 0{a^ In^) is the minimum order at which we 
can recognize the pattern of constraints. 
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B. Photon propagator 

Next we discuss the Schwinger-Dyson equation for the gauge boson. This equation has 
some different features from the fermion SDE. Now, the two fermion legs have to be treated 
equally. We can ensure this symmetry property by dividing the external momentum flow 
equally in the loop as shown in Fig. |4] : 

/ + q/2 




/- q/2 

FIG. 4: Unquenched Schwinger-Dyson equation for photon propagator. 
Using the Feynman rules, Fig. H] can be expressed as : 

which can be symbolically written as A~J(g) = ^(g) + 11^,^ (g), where II^j, is the photon 
self-energy and £+ = {£ + q/2) , £_ = {£- q/2). 

The definitions of the fermion and photon propagators ^ are given already in Sect. [Tll 

zSf{U) = iF{£^)/U , 

Eq. fl52|) must satisfy the photon Ward identity, q^A~l = q^q'^/C,, which is, of course, ful- 
filled by the bare propagator in Eq. fl52l) . Consequently, the loop graph of Fig. 4 must be 
transverse. Contracting Eq. (1321) with q^ and using the Ward- Green- Takahashi identity of 
Eq. ([1]), this transversality requires : 

Q,^'' = ^SfY, I dHi:T[Y{S,{U)-SF{£.)]] = 0. (53) 

^ Where appropriate, we denote the fermion and photon wavefunction renormahzation functions as F{p) or 
Fij)^) and G(p) or G{p^), respectively. Where we wish to emphasize that the quantities arc unrcnormal- 
ized, will be added to the list of arguments — with similar conventions for the rcnormalizcd quantities, 
for instance Ffi{p) and Guip)- 
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If dimensional regularization is used, then this integral is automatically zero. However with 
cut-off regularization, this is not the case. Then 11^,^ is not entirely transverse. To extract 



the correct component, we introduce the following tensor 46|, |47|] : 



Pf.u = ^ {^q^qu - q^g^^u) ■ (54) 

Projecting Eq. ( l52l) with P^^, allows us to remove the potentially quadratically divergent 
term in 4-dimensions, and project out the required ultraviolet logarithmically divergent 
terms. It is easy to check that this leaves the correct leading logarithms. We then have a 
scalar equation for the photon wavefunction renormalization : 

G{^) = 1+ N^'-^ F{nFih)P,uTr[T^^{i_,U,q)) n -(55) 

Recalling the definition of the vertex of Eqs. (171- [TOl) . we obtain : 



G{q\h?) 

xP,.| \l\ej%q') Tr(7^f+7'' 

+ Tr(r^ f-) j- (56) 

Moving to Euclidean space, we perform a Wick rotation. Substituting d'^i = 
2-71 i"^ (M! sin^ \1/ and the form of the transverse vertex from Eqs. (9, 10), and then taking 
the traces leads to : 

1 ^ _ aoNp f f df 



^(g^A^) QTi^q 



/•TT 

/ sin2^rf^F(£+)F(^_) 
Jo 

2\f{f,il,q^) (le^^l^- 3g2- 4£2 

I q^ 

+ riiiJj^^q') {2 {4f+ q') A^} 

+ ri{lJA^\q') {-8A2-3g2 [Af - q^)] 

+ r^{i-\U\q'){Ql-q{U'- q')] 

+ r3^(£_^£+^g^){8A2}l, (57) 
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where = ■ g)^ — i'^q^ and the photon Schwinger- Dyson equation picks out loop momen- 
tum regions where £4.^ ~ iJ^ ~ £^ ^ q^. This allows us to carry out the angular integrals in 
Eq. fl57p for the leading log terms after inserting Aj's and r^'s from Eqs. 1231) . This gives : 



1 + 



Qio Nf 
2,71 



3 

4'^ 



(58) 



where 



ry-^ ^ {S, + 62) r^'"^ - (53 + ^3) r^y^ + (56 + ^6) r/'"^ - (^8 + e,) r^^^ 



(59) 



This time the explicit longitudinal contribution comes from Ai; A2 does not contribute to 
the leading log's. Using Eq. 0411] and performing the radial integration, Eq. fISS]) yields : 



u=l 



+1 



A 

{U + 1) A2 



00 00 00 



-i E < E E E "o"«- 

u=0 m=0 n=0 r=0 



+m+2r+l ^^u+n+m+2r+l 



q^ 



-^7 

'-^ mrirr 



A2 (m + n + 2r + 11 



(60) 



where = (^2 + £^2) Sf^nrr - ('^3 + £^3) S^Lnrr + ("^e + e^e) '5^nrr ~ (<^8 + S^nrr 5 (61) 



2d 



(62) 



Evaluating the multiple sums using the symmetries and rearranging terms with respect to 
powers of ao yields : 



where 



1 



G(g2,A2) 



3 TT A"^ 

AT °° 2 



A.,,, 



A2 + 1^ 



b=0 

A'o = , 



(« + a-c+l) jrh)"^'^^'-'-'^ I '^<^(-'^)^^ ' 

^ ^ c=0 a=0 \ d=l ^ ' / 



(63) 
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with Rj defined by Eq. (150|) . Employing the expansion of the transverse vector coefficients 
introduced in Eqs. (12511271) . we can then write 1/G{q^) analogous to the fermion result for 
1/ FijP') of Eq. (I3T]) . after performing the many integrals : 



G(g2,A2) 



1 + 



Nf 
37 



olq In 



q_ 

A2 



t 



21 2 q 

CKn m -— ■ 



Oil ln=^^ 



Ai 3 - 



1000 



^2 

~3 2 1000 

A A ^1 '-'lOOO 



2000 



3 -^7 



0200 



0011 



^3 

4" 



4 



1 /li 9^ 

^ ^1 '^0011 



1 

2 ^1 '-'0200 



1 --y 

^ ^1 '-'2000 



+5, 



0300 



1011 



1^7 



1200 



3000 



(64) 



We have already made use of Eqs. (l30H32l) in above expression. Now the transverse vertex 
must have the right structure, i.e. the right coefficients A 
tion of the Schwinger-Dyson equations for 1/F(p^) and 1/G'(g^), Eqs. (jUl \5T[ 
multiplicatively renormalizable. 



'■mnrri ^rnnrri that the Solu- 



are 



IV. MULTIPLICATIVELY RENORMALIZABLE F{p'^) AND G{q^) 

A. The Photon Propagator 

We shall first look for the most general form of the multiplicatively renormalizable photon 
wavefunction renormalization. In order to do so, the renormalized Gr can be written in the 
following form by using Eq. f|T6l) : 

G^{q\^^') = ^3-l(/.^A2)G(g^A2). (65) 
23 



We define the most general leading logarithmic expansion of the unrenormalized photon 
wavefunction renormalization by : 

oo 2 



A2 

u=0 



2 ^2 ^2 



= 1 + aoSi ln-|^+ alB2 \n^j^ + alB^ In' + O(a^) . (66) 

We impose the renormalization condition that GR{q^ = /i^) = 1. The coefficients Bi {i > 2) 
are then constrained by multiplicative renormalizability, i.e., B2 = B^, Bn = (-Bi)" so that 
the renormalized photon wavefunction renormalization can be written as : 

00 2 

u=0 

2 2 2 

= 1 + a« 5i In ^ + al B^ In' ^ + al B^ In' ^ + C(a^) . (67) 

^[X ^[X ^t/ 

Then, as we shall use later, the inverse of G and Gr are : 



1 q2 1 q2 

1 - ocrBi In— , „ = 1 - ocqBi In— 



(68) 

where, as is well-known, 1/G(g') in QED only has a leading logarithm at one loop order, 
just like l/a{q^) and being related to this physical quantity is independent of the gauge. 



B. The Fermion Propagator 



Analogously to the previous section, we deal with the fermion wavefunction renormalization. 
We similarly define the general leading logarithmic expansion of the unrenormalized F as : 



A^ 
2 

P ,^24 i„2 V" , 3 . ,3 P" , 



A 

u=0 

,2 J2 J2 



1 + «o In ^ + A2 In^ ^ + M In' ^ + C?(«o) • (69) 
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Since not only the coupling, but the gauge parameter have to be renormalized, we need to 
make the dependence of the Au on ^ explicit. As gauge dependence in the coefficients arises 
from photon propagators, any Au cannot have a higher power of ^ than C,^. Consequently, 
FijP', h?) can be written as : 

F{p\K^) = l + ao(aie + &i) In ^ 

2 

2 

+ ao'(«3 e + he + czi + ds) In' |^ + 0{a^,) , (70) 

where a,, 6j, Cj, di are constants related to the A^ by comparing Eqs. (169|) and (170|) . Recalling 
Eqs. ( HTt dn]), ^0 = Z-s^R , ao = Z^^ we note that 

and = Z2"^(/iVA')i^(p',A'), (71) 

with the renormalization condition for the fermion wavefunction renormalization 
FjiijP' = /i^) = 1. Eq. (1701) can then be inserted in this equation to give 



/W^) = 1 + "/? (oi^R + h) In — 

+ «R(a2 + 62 + C2) In^ — 

2 

+ alias + &3 e« + C3 + 4) In' ^ + 0{aj,) , (72) 

Multiplicative renormalizability requires that the inverse unrenormalized fermion wavefunc- 
tion renormalization must have the following form keeping only the leading logarithms : 
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F(p2,A2) 



1 + ttoln— [-ai^ - 61] 



+ c^l In^ P 



2l^2P_ 

A2 



4i^4P 



2 r 
2 

A2 



6 2 ^ 



2 ^ ^ ^ 6 3 



24 6 4 ^ ^ 

4 V 6 ^ 6 ^ ^ 



h\B, 



(73) 



The renormalized form of 1/F can be found by replacing — > — > E,r and A — *■ /i 
in the above expression. 



V. MR CONSTRAINTS ON THE VERTEX 

In Sec. mil we have shown exactly how the full vertex contributes in the fermion and boson 
SDEs. In principle, truncation of the Schwinger-Dyson equations for the fermion and boson 
propagators requires knowledge of the complete structure of the vertex, all 12 independent 
components or here in massless QED all six. While two are fixed by the Ward-Green- 
Takahashi identity in terms of the fermion propagator, the four transverse components 
appear to embody information about all the higher point Green's functions. Knowledge we 
do not have, unless we solve the theory completely. However, two simplifications have already 
occurred. Firstly, the massless fermion and boson self energies involve just two projections 
of the six independent vertex vectors, so we do not need to know their complete spin and 
Lorentz structure. This is helpful, since even at O(a;o) in perturbation theory, this is of 
daunting complexity [41|. The second simplification is that multiplicative renormalizability 
(MR) is closely related to the ultraviolet behaviour of the loops in Figs. [3|, |H There not only 
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is the structure of the vertex simpler, but importantly for the present study the 2 graphs 
explore the vertex in distinct kinematic regimes. For the fermion self-energy, the internal 
fermion momentum k and boson momentum q are very much larger than the external fermion 
momentum p, i.e. k"^ ~ >> p^. In contrast, for the boson self-energy, it is the fermion 
momenta that are both large, i.e. /c^ ~ >> g^. We shall see that this distinction plays a 
powerful role in our analysis. 

First, in this section we combine the results of the previous two sections to find the con- 
straints on the fermion-photon vertex imposed by multiplicative renormalizability. 



A. MR constraints via fermion Schwinger-Dyson equation 



In this and the next section, we apply the above strategy first to the fermion wavefunction 
renormalization in full massless QED. To do this, we start by comparing order- by-order the 
results fixed by multiplicatively renormalizable F, Eq. fl73|) . with those found by solving the 
Schwinger-Dyson equation, Eq. (ISTl) . These comparisons will give what we refer to as the 
fermion conditions, labelled by FCl, FC2, etc.. 



cxq In—— comparison: 



-Ai = -{ai^ + bi) = - ^ 

4 71 

ai = , bi= 



(74) 



In this first order comparison MR fixes the value of ai and bi and by that all leading order 
terms in 1/F or F, then Eq. (173|) requires 



FCl 



A^= ^,A,= ^ 
47r' 2 



(75) 



21 2 P 

a.Q In comparison: 



2 ^ ^ 2 47r 



^ 3\ , 3-i-/ 3-/ 



1000 



(76) 
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Making use of Eqs. (17^ 1751) and keeping in mind that Aiqqq and iS^goo '^^^ be at most 
proportional to ^ or Np from Eqs. fl2Sl we immediately see that the ,^^-term on both 
sides automatically matches and for the ^-term we must have : 



FC2 : 



2 ~ -^1000 '-'1000 



(77) 



3i SP 

ckq In comparison: 



3! 



1 

4:71 



3 



1 Ao 



Ai Bi 



+ \{A, + B,) [a 



1000 '-'lOOO 



A-^ 

■^2000 



1100 



-4^ 



2000 



0200 



— o, 



0011 



(78) 



The leading terms in C, in Eq. flTHj) {i.e. 0{^^)) automatically match on the left and right 
hand sides. Imposing Eq. (1771) . the 0(^^) terms require the transverse part to be fixed so 
that 



FC3 



4 



-A. 



2000 



1100 



4^ 



2000 



4^' 



0200 



0011 



(79) 



As one can see the Bi term in Eq. (!78l) disappears from the above expression and this must 
repeat itself in every order, i.e. in leading order terms the photon contribution will be 
cancelled out by the transverse vertex. 
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4i 4P 

a.Q In comparison: 



24^ 4! Att 



A1A2 ( - + - e 



4 ^ ^ 16 



— Ai B2 
16 



+ -A^ 



1011 



+ {Ai + Bi)(-A, 



^2000 
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S 
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- Al 5l (a(qqq - S 



'-A' 

o '^1100 
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3 9/ _ i 9/ 
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is' 



1200 



3000 



(80) 



Once again in above expression the leading terms in ^ (i.e. 0(^^)) terms match on both 
sides. After substituting the FC2 and FC3 conditions in Eq. flHOj) . we have the following 
combined constraints on and Np terms : 



FC4 : 



16 


(Ai + B,) 
48 


^A^ - 

■^1100 


'5 0011 + 


2000 ^ 


'^0200 1 






+ + 

' g -^1011 ^ 


^a' 

-^1200 


+ -A^ 

' ^ -^2100 


3 

+ i 


A^ 

•^3000 






+ ^^oni + 


^ '-'0300 


12 


1 

+ 6 


'-'1200 


2 '-'3000 • 



(81) 



The Schwinger- Dyson Equation for the fermion propagator involves corrections from pho- 
ton emission and absorption as displayed in Fig. HI This requires the fermion renormal- 
ization function to depend on the photon renormalization function, which in turn de- 
pends on the number of fermions Np. Therefore in general | ^noo' '^ooii' '^/ooo' '^020o| 
{Moii^M20o^Mioo^Mooo^^oiu^^Lo^^Li^^!20o^Lo} ^rms in Eq. §3) can be propor- 
tional to {^"^ or N| or C,^f) and {C,^ or C,^Np or or Np) respectively. Remarkably, the 
matching required by multiplicatively renormalizability of these renormalization functions is 
automatically satisfied if the transverse fermion-boson vertex is independent of the photon 
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renormalization function at leading logarithmic order. Therefore | A(iqq, ^qq^^^, iS/qoo? '^0200 | 
and |Ami'A^2oo'^2wo'A^oo'^oni> ^0^300' ^K)ii'^4oo ^rms would be proportional 



to only and C,^ terms respectively. This will clearly constrain the non-perturbative forms of 
the transverse vertex that we wish to determine. In other words constraint FCA of Eq. (IHTI) . 
will divide into two separate conditions for ^"^Np and C,^ comparisons : 

FCAl : = AiiQQ — Sqqh + 1S2000 + '^0200 ; 

FC42 ■ ^ - -A^ + —A^ + -A^ + -A^ 
^ '-^^^ • 16 ~ 6 12 ^200 ' ^ -^2100 ' ^ -^3000 

12 4 0300 -|^2 1011 ' Q '^1200 2 3000 • \°'^J 

The idea is then to find a non-perturbative structure for the transverse pieces that delivers 
such relations. This we do in the next section. However, first we determine the conditions 
imposed by multiplicative renormalizability for the photon wavefunction renormalization. 

B. MR constraints via Photon Schwinger-Dyson equation 



We now repeat the previous steps for the photon wavefunction renormalization. Compar- 
ison takes place between Eq. flMl) and Eq. fpSl) order-by-order for 1/G. Obviously, this 
time instead of looking at the terms depending on the gauge parameter ^, we compare the 
dependence on Np, the number of flavours hidden in the Bi terms. These give what we refer 

to as the photon conditions labelled PCI, PC2, etc.. Then : 
^2 



1 P 

ao In-— comparison: 



PCI 



Nf 
37 



3 7r 



(83) 



First order comparison defines the value of Bi in terms of Nf and as given in Eq. fl67j) fixes 
all the higher order terms. 
,2 



cKq In^-^ comparison: 



PC2 



-A. - 



(84) 



As we see above the second order comparison imposes this condition on the symmetric part 
of the transverse vertex. 
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A2 



comparison: 

6 



^-^7 _|_l"c')' M-'^'^ 
2 1000 + ^ '-'2000 ^ '-'0200 + ^ '-'ooli 



Substituting Eq. ( 18^ in above condition yields : 
PCS : 
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1 ^ 
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3^ 



0200 



2 '-'2000 ' 



where every term is proportional to ^"^ 
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-^1 ^1 "7=* 7 ^1 "7=» 7 ^1 -pr 7 ^1 7 1 7 

= ^ '^0011 ^ '^0200 + ^ '^1000 + ^ '^2000 + g '^3000 
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1200 



(87) 



Making use of Eqs. 



PC4 : 



the above expression becomes : 



24 
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3000 5 



where every term is proportional to C,^. So far we have expressed the general multiplicative 
renormalizability constraints on the 3-point vertex function in terms of the constants A2,nrr 



and S2,nrr up to 0{a^). 



C. Generalized Fermion and Photon MR constraints 



Let us first look at the general picture. Firstly, ai and bi being fixed by Eq. (!74|) . allows the 
expansion coefficients Au in Eq. (l69l) to be fixed in all orders : 



f 42 43 

47r' 2! ' 3! 



4 ^ 

A - ^ 

Ul 



(89) 



and in turn the infinite leading log series of F{p'^) in Eq. ( l69l) can be summed up as a power 
series : 



u=0 



Au „2 
A2 



2 \ aA\ 



P_ 
A2 



(90) 
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This is the non-perturbative expression for the unquenched (full) fermion wavefunction 
renormalization. Moreover, it has exactly the same form as in the quenched QED 

49I. Secondly, the relation between the photon coefficients 



1 Cll VJ l lll d llZjCl 

Q, Q, 
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are also found through PCI 



Nf 
37 



Nf 
37 



hence the infinite series of G{p^,A'^), Eq. fl66|) can also be summed up as : 

.2 



u=0 



aoBi Ing2/A2 ' 



(91) 



(92) 



1. Generalized MR constraints from fermion SDE 



Making use of Eqs. (1551 - [92]) . we can then rewrite the inverse fermion wavefunction renor- 
malization calculated from SDE, Eq. (1471) as : 

1 f 1 



F(p2,A2) 



1 - <^ - 



F(p2,A2) 



+ 1 



it=i 



2 ^ a! ^ M + 1 

a=l 

1 " ^ 4'' 4° 1 

_ \ \ ^_^ \b tJ±tJ± r>u-b-a 



b=l a=0 



M = l 



(93) 



and as a consequence of equating the multiplicatively renormalized F, Eq. (1901) to Eq. (!93|l 
we can extract the generalized MR constraints to all orders, which, of course, reproduces 
FCl to FCA : 



00 2 
^ ao"+^ ln"+^ 1^ (i/„ + F„) 

«=i 

00 2 f 1 " 1 

M=l L a=l ^ 
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u u—b 



'6— a 



b=l a=0 



+ 11 



(94) 
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where in Eqs. ( HHl H9l) for the Hu and one can now substitute for the An-, Bn from 
Eqs. dMlEIl). 

2. Generalized MR constraints from photon SDE 

Making use of Eqs. fl89ti92p . we repeat the above procedure for photons, which is analogous to 
the fermion case above, in order to rewrite the inverse photon wavefunction renormahzation, 
Eq. ([63D : 

1 Nf g2 

-S|:«rin-|;{(^.^A,}, (95) 

where in the expression for of Eq. fl63|) we can substitute the conditions for An from 
Eq. (!89|) . The generahzed MR photon constraints can then be written as : 



(96) 

automatically satisfying PCl-PCA. 

D. Non-perturbative Fermion and Photon MR constraints 

1. Non-perturbative MR constraints on transverse vertex from photon SDE 

To understand the above conditions in full generality {i.e. beyond their expansion in leading 
logarithms) we first turn our attention to photon equation. Starting from Eq. (!68l) for 
multiplicatively renormalizable 1/G{q^, A^), we see that multiplicative renormalizability for 
leading logs is given by just the 0{ao) term, Bi, which is gauge independent. Equating 
the photon Schwinger- Dyson equation at the leading logarithmic order, Eq. ( 158|) . with the 
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multiplicatively renormalizable G{q^), Eq. (168|) : 
1 ^ a^N, 



_|_ _j7sym 

(£) 4 ^ 



We observe that the Ai term of the Ball-Chiu longitudinal vertex generates this. However, 
importantly for the present purpose this is part of a whole series: 



1 



2 6 24 120 720 



5040 



40320 



(98) 



where X = oqAi In—. Beyond 0{ao), this series {i.e. terms inside the square bracket) 
has to be cancelled exactly by the contribution from the vertex components. Since the 
A2-term in the Ball-Chiu longitudinal component only contributes at non-leading order, it 
is the symmetric part of the transverse vertex, r^^™, with its implicit gauge dependence 
that has to provide this cancellation. PC2 to PCA in Eqs. ( IMHHSl) give the conditions for 
this cancellation to be achieved at 0(0^), O(tto) and O(ao) and the general condition in 
Eqs. ( l96l 1971) for all orders. To go further, we note that multiplicative renormalizability of 
the photon Schwinger-Dyson equation, Eq. fl971) . picks out loop momentum regions where 
£^ 2^ £^ ~ £2 The second term in Eq. (197|) must give the following result : 



Air 
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Nf , q' 
3 vr 



Ix + -x2 + lx3 

2 6 24 
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— X' + ^X^ + 
720 5040 



40320 



x' + o(«^ 



(99) 



This surely determines the structure of the r^^™'s for this to happen. The dependence 
on the fermion wavefunction renormalization must be more complicated than 1/F times a 
kinematic factor. It must be proportional to a function of a function of F's so let us write : 

1 



— sj/m 

7 



F{q) 



h{Y) 



(100) 



In keeping with the ethos of this work, we assume that Y is determined by the fermion 
wavefunction renormalization. Since the renormalization of the Tj's is replicated wholly by 
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the factor of l/F, Y must be renormalization independent. As an example let us choose it 
to be: 



Y = ^-l (101) 



where the factor of —1 ensures that the leading logarithm expansion of Y begins at 0{aQ In) 
as required by Eq. (l99l) . Can we find what function h{Y) is to satisfy Eqs. (1971 1991) ? Let us 
assume we can expand h{Y) as a power series in F, and in turn expand this in leading logs 
of momenta. Then to produce the cancellation required, we deduce : 

MY) = Y + -Y^ - - Y^ + — Y^ - —Y'^ + — Y"" ~ —Y' + 1^8 + 0(^9), 
^ ^ 2 6 12 20 30 42 56 ^ ^ ' 

(102) 

We recognise this as 

-YY 



^ nin + 1) 

n=l ^ ' 

l + F)ln(l + r), (103) 



substituting Y from Eq. fllOll) . h{Y) becomes : 



Mn-^m^ . (104) 



Since a form like r^^™ ~ l/F{l'^) \n{F{q'^)/ F{f)) in Eq. ffTOOll is the k ^ p = i limit of 
the evolving structure, this naturally generalises to the k ^ p configuration as : 

While the form in the photon limit is determined, the structure in general momentum 
configurations is not unique and there are several possibilities differing only beyond leading 
logarithmic order. Three of these are: 



2 \F{k'^) F{p^)J 2 \F{k'^) F{p 
5(^)^1 ' 1. ^C')^ 



2 (F(A;2)F(p2))i/2 F{k^)F{p^) 
1/1 , 1 \ ^ F(g2^2 



5(3) = t ^_ + ^_ In _LM2 (106) 

4 \F{k^) F{p^) J F{k^)F{p^) ' ^ ' 



all of which give the same h{Y) of Eq. (11041) in the photon limit of A;^ ~ ^ g 
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2. Non-perturbative MR constraints on transverse vertex from fermion SDE 



Similarly, for the multiplicatively renormalizable 1/F(g^, A^), the result at leading logarith- 
mic order is given by the leading ^ dependent piece, as required by the Landau-Khalatnikov- 
Fradkin transformation 43|, . This leading term is provided by the first term in the inte- 



grals of Eqs. (l43l l45l l47l) . Let us recall Eq. (l43l) and in this equation we perform both the 
radial and angular integration for the first term, but only the angular integration for the 
second term, then we find : 



1 



1 + 




1 



.(107) 



Imposing the MR fermion condition, Eq. (I90l) . on this expression yields the following con- 
straint on the transverse vertex : 



3 tto 



1 



1 



2 yF^k"^) F{p^ 



0. 



(108) 



This cancellation involves both the longitudinal and transverse pieces together. At leading 
logarithmic order the longitudinal contribution comes from just the A2-term in the Ball-Chiu 
vertex. 

While antisymmetric forms do not contribute to the leading logarithmic behaviour of the 
photon Schwinger-Dyson equation, this is not the case for the fermion equation. Indeed, 
here the distinction between symmetric and antisymmetric disappears when /c^ ~ ^ p^. 
Thus, a seemingly symmetric form like 



In 



F{k^) F{p^ 



In 



F{k^ 



aoA^ln— + 0{al), (109) 



is antisymmetric in k and p . Such a form contributes equally to the antisymmetric terms 
like 



1 



k^ 

-aoAi In — + 0{al) 



(110) 



F(P) F(p2) 

The 0{al), O^a^) and O^a^) conditions of Eqs. (I84ti88p . which embody the gauge indepen- 
dence of the photon wavefunction renormalization and the known gauge dependence of the 
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fermion function, require the transverse vertex to deliver a very particular gauge dependence 
itself. Our aim is to reproduce this by constructing the non-perturbative transverse vertex 
from the fermion wavefunction renormalization. This means from Eq. (11081) that 

:=anti , :=sym _ _1 ( /iiiN 

Hence this expression tells us that the total transverse vertex, i.e. combination of anti- 
symmetric and symmetric parts, must be proportional to antisymmetric form in the limit 
/c^ ^ p^. These considerations suggest particular antisymmetric and symmetric ver- 

tex forms. In Table I, we give the specific coefficients Amnrr and Smnrr at C(ao) for these 
examples. 



VI. APPLICATION 

The next step is to make use of all the examples in Table [T] as inputs to the multiplicative 
renormalizability constraints. In order to satisfy these, we have a set of equations to solve. 
As a first step the coefficient functions, r^'s, can in general be written as a sum of different 
non-perturbative forms of F and G using the above examples. Hence, an antisymmetric and 
symmetric combination of F and G in rf"** and r/^"^ respectively become : 

^sym ^ IJH) ^(1) j ^ IJi2) ^(2) j + . . . + (^/(n) ^(n.) ^ ^^^^^ 

where A^"'^ and 5'-"^ refer to the relevant expressions in the left hand column of Table [B 
In general the number of constants needed to solve these equations is proportional to the 
number, n, of various combinations of the F and G. These combinations will appear in the 
ansatz for the non-perturbative transverse vertex. We then try to solve these equations by 
choosing a minimal number of combinations, in order to find the simplest possible vertex 
ansatz. 

From Eqs. (l46l [611) . we see that the coefficients 7^, ej, defined in Eqs. (!25|) appear in 
the fermion and photon conditions in rather specific combinations. To make this explicit 
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1 1 

F{k) F{p) 
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5(1) 


1 ( 1 
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c /iS c _ "^1 c _ "^1 

03000 —-777^1) 02100— — r, 01011— -777 
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5o30o =5i2oo= 5oiii=0 
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^F{k) F{p))'''F{k)F{p) 


01000 ~~~2' '-'0100 — ^ 

A^ A\ A\ 

^2000 =— 5iioo= —-^1 5ooii= — 5o2oo=0 
4 2 4 

^1 c ^1 c ^1 

03000——-^) 02100— -^1 01011— - — 
5o30o =5i2oo= 5oiii=0 



TABLE I: Antisymmetric combinations of F and G. 
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and simplify the notation, it is useful to define : 

Pf = {P2 + + - Ps) , 7/ = (-73 + 76) , 

6f = (62 + S3 + 6g - 6s) , e/ = (-eg -eg) , 

^7 = (62 - S3 + Sq - 5s) , = (62 - £3 + ee - es) • (113) 

Recall that antisymmetric forms for the Tj's do not contribute to the photon renormalization 
at leading logarithmic order, and so we have no corresponding combinations of and 7^. 



A. Fermion constraints 

We now wish to write down the fermion constraints FCl — FC4, Eqs. (1771 ISTl) . which we 
obtained in the previous section for the specific choices for r^"** and r^^™", namely A^^"^ 
as the antisymmetric form of the transverse vertex and S^^^ as the symmetric one in the 
Table [T]: 

^(1) = L \ 5« = i + ^] In ^ + ^] (114) 

V^('^') Hp^)) 2 \F{k^) F{p^) J 2 \F{k^) F{p^) J ^ ' 

After recalling the definition of A^^nrr 5^^^^ from Eq. (jlH]) and reading off the specific 
values oi Amnrr and Smnrr's from the TablelH the MR constraint FC2, Eq. fl77|l . which comes 
from In^ p^/A^ order comparison together with Eq. flllSp gives the following condition : 



-{/3f + lf)+ ^{Sf + Sf) 



(115) 



The In^ p^/A^ order constraint FC3, Eq. (1791) . splits the combined /?/, 7/, (5/, form of 
previous constraint into two separate ones : 



(Sf + Sf) = 


0, 


iPf + ^f) = 


1 

~2 ■ 



(116) 



The In^p^/A^ order constraint FC4, Eq. 
but again yields Eq. f lll6p . 



(IHTj) does not give further new information, 
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B. Photon constraints 



We repeat this procedure procedure for the photon constraints PC2 — PCA for the same 
choices of ^^^^ and 5^^) in Table [11 All the MR constraints PC2 to PCA, Eqs. dHl [861), 
which follow from Iv? to comparisons give the same condition and that is : 



{5. 



'7 ' ^7/ 



4 
3 



;ii7) 



Since this condition repeats itself at every order, this means we have the exact solutions! 
There are 14 constants to be fixed, and Eqs. (11161 (I117|) can only fix three of them in terms 
of the others, for instance : 



1^2 = -- - /^s - /^e + /^s + 73 - 76 , 

02 = - Oe + Os - — + - Eq + — 

A = - + ^ - ^ 

^ ~ 3 2 2 ■ 



(118) 



Substituting these constants into the Eq. fl2^ we can write the non-perturbative coefficient 
functions r,'s as : 
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(119) 
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with T"""**' having specific forms such as those determined in Sec. V, Eqs. (11061 II 111) , 
examples of which are given in Table [H 

Multiphcative renormahzabihty relates the coefficients at order (ag In)" to that at n = 
1. This lowest leading logarithm coefficient is fixed by the longitudinal component of the 
fermion-boson vertex. Transverse components only enter at = 2. Remarkably, once 
the MR conditions at this first non-trivial order are satisfied, the conditions at all orders in 
leading logarithms for both the fermion and photon Schwinger-Dyson equations are fulfilled. 

As far as the leading terms are concerned, the above constraints ensure that both fermion 
and photon propagators are multiplicatively renormalizable in massless unquenched QED4. 
These constraints impose conditions on the transverse part of the vertex. The 3-point vertex 
calculated at O(ao) and the coefficient constants, r^'s, at one loop order 4l|] will be very 
helpful in fixing some of these constants. 



VII. PERTURBATION THEORY 



The vertex coefficients r^'s were calculated exactly in O(a;o) for the massive fermions in 



a general covariant gauge 4l|] and for our purpose their massless limits are given in Ap- 
pendix El 



We observe in Eq. flAltlA4|) that all the four Tj's (z = 2,3,6,8), contains four different 



structures in general. The first one is the Jo dependent part, which contains Spence functions 
(or Dilogarithms) of momenta p'^,k^,q'^ in Eq. ( 1A5I) . The second part is proportional to 
In/c^/p^ which is the perturbative expansion of the asymmetric combination of F and G 
in first order, and the third part is proportional to lng^/(/i;^p^), which is the perturbative 
expansion of the symmetric combination of F and G, and the final one is the kinematical 
term dependent on k^,p'^,q'^. 

In order to fix some of the individual constants 7j, Si, e^'s appearing in Eqs. (11161 1117]) 
we need to make a comparison between perturbative transverse vertex coefficients rf^*"* of 
Eq. (IAH - IA4|) . and the non-perturbative ones we used in fermion and photon SDE, 'j-^°^-p^^*- 
of Eqs. ( l26l 1271) in the previous sections. However this comparison has to be made in a 
particular way in order to be meaningful. There are two points to be considered. The first is 
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how these coefficients behave inside the fermion and photon SDEs, since these equations 
project out different parts of the vertex. Recall, that with this in mind we started with a 
simplified ansatz for the explicit kinematic factors in the ^-J^"""?"^^*-^ ([53]), and assumed 
their denominators did not depend on k ■ p. We therefore need to take the corresponding 
limits of both pure perturbative rf^^^ 's, Eq. (lAlj - IA4p . and the -;-^"°"-P^'*-'g^ Eqg^ [27]) 
which we inserted into SDK. While for the fermion SDE the relevant limit would be where 
either of the fermion momenta are large, e.g. k"^ c::^ q"^ ^ k -p ^ p"^, for the photon SDE the 
relevant one is where the both internal fermion momenta are same and much greater than 
the photon momentum, e.g. /c^ ~ ^ q^. 

The second point is that the real functions depend on the angle between momenta k and 
p. This means that when we obtained MR constraints, Eqs. (111511117]) . on the vertex, i.e. on 
Ti functions, their angular dependences were already integrated out. These angular averaged 
functions we call effective r^'s jSO]]. It is these that we have to compare with perturbation 
theory. 



A. ^ p^: The Fermion Limit 

Let us take the fermion limit of the perturbative rf^^* 's, Eqs. flAlj - IA4I) in Euclidean space. 
In order to do this, Jq of Eqs. (]A5l IB1|) has to be expanded up to O^l/k"^) to ensure we keep 
all the terms of the required order. As shown in Appendix [B] these results are : 

/ Expert. / 2 , 2 2n "O ^ , J 4 k ■ p U p'^ ^ 

(^E-.pert. , 2 ,2 2^ _ ^0 \r, — { - M ^lll. + A ^ \ 
/ E\pert. / 2 I 2 2\ Q^O ^ , f 1 1 k ■ p 1 p^ 

(^6 Yneai (P ,k ,q) = —— In^J 



87rfc2 p2 i_ 3 3 A;2 5 ^2 
= 0. (120) 

In this limit one observes that both Jq and \n. {q^ / k'^ p^) behave like \n. [k"^ / p^) . Therefore 
all four coefficient functions become proportional to In {k"^ /p^) signaling that the structure 
of non-perturbative transverse vertex consists of purely asymmetric combination of F or G. 
Next we expand the non-perturbative ^^"°"-P'=^*-'g^ ([25]) . using Eqs. (!28l [29]) at the order 
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As we mentioned earlier, during the process of finding MR constraints in Eq. ( 1116^ from the 
fermion SDE we performed both radial and angular integrations therefore these constraints 
on the vertex are for the r^'s whose angular dependence has been integrated out, viz. they 
are the effective ^-"'"^"P'^^^ 's. To make consistent comparison between the Eqs. (112011121]) . we 
must integrate out the angular dependence of both rf^*^*' and ■jr^°"--P^^^- ^ xhe details of this 
procedure can be found in Appendix O Following this, the effective coefficient functions 



can be found from t^J^i's in Eq. (11201) : 
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We repeat the same procedure for the first order expansion of the non-perturbative coeffi- 
cients 7"]Jear^'^'^* ^'i- P^ip to give : 

2 / P\ 



1 . / , k 



(^^^)non-,eH. ^^2^ ^2) ^ ^ ^ ("^o ^00 In ^ 

+ p(^6 + e6) Uo^i'oooln^ J 

(ts^):;;-"'^' (p"^ = ^ /^^ (-«o^?ooo in ^ 

+ ^ '^s (ao5,«oooln^) ■ (123) 

The constants /5j's, (5j's, 7j's and £j's appearing in Eq. (11231) are the ones which must satisfy 
the MR constraints, Eqs. (11161 11171) . Let us check we have obtained the correct result in 
three key situations. 

First we compare Eq. (11221) with Eq. (I123p to read off the constraints on ^{ooo ^nd 5]*ooo 



i = 2,3,6,8 



-4,1000 ~ '^2 ^1000 — 



{^3 - Is) ^?ooo - (<^3 - es) 51*000 = ' 

iPe + 76) -4iooo - + ^e) ^fooo = ^ > 

PsAUo-SsSUo = 0. (124) 

la) General Case at O{ao) : Recall the definition of ^looo ^nd iS^ooo' ^'i- P6|^ - in or- 
der to form the FC2 constraint in Eq. (I77|) using above expressions by adding them up 
appropriately : 

a(ooo - 5/000 = - 1^ + i) ^1 = T • (125) 
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lb) For the special vertex (A'^^^ and 5^^^) at C(ao) : Making use of Table [T] we can 
read off the value of A\qqq and SIqqq and insert them into Eg. (11 241) to see whether we 
can satisfy the fermion MR constraint of Eq. (IllSp by using Eq. (11130 : 

m + + Ps) + (-73 + 76)] (-^1) - [(^2 + S3 + 6e- Ss) + (-£3 + e,)] {-^) 



-if3f + lf) + ^{Sf + ef) 
3-12 + 4^- 



A 
2 



(126) 



As we see all effective T^f/s, Eq. (11221) add up to Ai/2, as required. 

2) Non-perturbative check : If we trace back the MR constraint in fermion SDE equa- 



tion, Eq. (11081) . we have already observed that the ^ dependent part will give the right 
equality and the rest must be zero to give the fermion MR condition. Hence this MR con- 
straint for the effective r^'s after the angular and before the radial integration was performed 
can be written as : 



3ao 



57r 
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2 \F{k'^) F{p'^ 



0, 



(127) 



where 
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Making use of Tg^j's in Eq. (I122p to form Eq. (I128p gives 



ao In 



p2 



2 1 1 
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Since Eq. f lM cancels out in Eq. ffTm?]) . Eq. ffT^ is satisfied. 



(129) 



(130) 
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B. A:^ ~ » : The Photon Limit 



Let us turn our attention now to the photon hmit of the perturbative rf^'^*''s, Eqs. f lAll - IA4p 
in Euchdean space. The technical details of this limit can be found in Appendix IB 21 and 
then we have : 

irifrlaii^^l') = +0{al), (131) 

since in this limit In {k'^/p'^) approaches 1 and In (g^/ {k"^ p^)) approaches In (g^/^^). Therefore 
all four coefficient functions become proportional to In (g^/ 1'^). This signals that the structure 
of the non-perturbative transverse vertex consists of purely symmetric combination of F or 
G. We expand the non-perturbative ^^"""-p^'^^ 's^ Eq ([^ . using Eqs. at the order 

O(ao) : 



]^ 02 



I) , 



t q 

{riTrTar^H^^q') =-]-, (58 + ^8) ao5«ooo In ^ + • (132) 

Comparing Eq. (11311) and Eq. (11321) one can read off the symmetric coefficients as : 

(<^2 + 82) '^1000 = ; 

^1 

(^3 + £^3) '^1000 = ^ ) 

(^6 + e^e) ^51000 = 1 

{5, + e,) 5fooo = . (133) 
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Analogously to the fermion case, we now perform similar checks for the photon constraints 
in the same three situations : 

la) General Case at Oi^ao) : 

Recalling Eq. (IM]) let us check whether the photon MR constraint PC2, Eq. ([HID, ^i<^o) 
is satisfied by Eq. (11331) after adding them appropriately : 

(^2 + 62) Slo,o - (^3 + ^3) 5fooo + (^6 + e,) 5i^ooo - ('^8 + es) 5fooo = - (^-^^ ^ , 

t.e. Si,, = ^A^. (134) 

lb) For the special vertex (^(^^ and 5^^^) at O{ao) : 

We also check if the photon MR constraint, Eq. (11171 ) at O{ao) is satisfied for this special 
choice of the vertex : 

((52 + ^2) -{S-i + e^) +{5e + ee) -(^g + ^s)) (^) = 



. 2 y 3 
4 

i.e. + = — . (135) 

3 

As we can see from both results, Eqs. (11341 1135]) . the effective t^j-j-'s satisfy the photon MR 
constraint. 

2) Non-perturbative check : 

Recalling Eq. (pTIl and after extracting the non-perturbative MR constraints, we can usefully 
rewrite this as : 



E\non—pert. / n2 2 



{i\q')}=0, (136) 



3 7r J„2 £2 ] ^ ^ 2 ^ ^ ^ ^eff 

where 

J2(rf):]]-^-'- if, q') = f {T2%ff + (riUf + {ri^f . (137) 

At O{ao) 

/f ^ [F{f) -^]=-Y "0 J-2 + ^(«o) • (138) 
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Making use of Eq. fll3ip to form Eq. fll37p we obtain : 



2 



ao In' 1^ + 0{al) 



(139) 



We see Eq. ffT^ cancels Eq. ffT^ and so Eq. 



is satisfied. 



C. Individual coefficients 

With guidance from perturbation theory, we can now find further relations between the 
constants, Eq. fll24p and Eq. fll33p . These eight equations fix eight of the 14 unknown 
constants (^2, ^s, ^e, ^s, )• In general these are : 





2 t2'-'i000 ' 


>^1000 ~ 


-Ai - £2 '5^000 ' 


^3*5^000 = 


^1 ^ .3 
2 fc3^1000 ' 


{03 - 73) ^?000 = 




^6'5|'qoo = 


^et-'iooo 5 


iPe + 76) ^?ooo = 


1 ; 




1 8 ^1000 > 


/Ss-^iooo ~ 


' 8 
~^8 '5^000 • 



(140) 

For the specific choices of antisymmetric, A^^^ and symmetric S^^^ transverse vertex forms 
given in Table [H Eq. fll40p becomes : 

2 1 

-3 + ^3-2^8, 

^+73-^3, 
1 

-^-76, 

1 

(141) 



^2 = 


2 

-- + 2^3 -eg , 


P2 


^3 = 


2 

3-^3 




^6 = 






58 = 


—£s , 


Ps 


£2 = 


-2 53 + eg- 





As we can see the unknown constraints in r^'s, Eqs. (125 | I119p . have now been fixed to match 
with perturbation theory. If we insert these constants in Eqs. (l25l I119p . we can write the 
coefficient functions, r^'s, in Euclidean space to obtain our final non-perturbative result : 
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where 



2 

1 

1 

1 

^ ~(P+p2) 

{k^ + p^y 
1 

(/i;2 — p2^ 

1 

~(A;2+p2) 



and Tj 
OR 



2 , £8 
o + ^3 - T7- 
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2A; -p 

fc2 + p2 



. anti 



{k + pf 



- £3 ) + 73 



2 g 
77-^3 



A;2 + p2 

1 g2 

- 76 



k"^ + p2 



sym 
^2 ) 



. anti 



4 '"^+^2 
,2 



anti 



-^6 



9 



/c2 + p2 



sym 



1 2A;-p 



anti 



-£8 



k"^ + p2 



. sj/m 



F(A;2) F(p2) 



1 / 1 



1 \ 



4 '^F(A;2) + F(p2) ) \F{k^)F{p^) 
I f I 1 



F{q 



2\2 



In 



1 /F(g2) F(g2)\ 

2 Vi^(P) F(p2) j 



(142) 



2 Vi^(A;2) F(j92) 

The fermion and photon SDE's at leading log order do not fix the constants 7^, e^, Eq. (11421) . 
As the simplest example for later exploration we choose 7^ = = in the above expressions 
and insert the second form of r/^'" in Eq. (I142p . we then have : 
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1 



+ 



- p4) 
1 

1 

1 



1 



1 



5 

12 

1 

3 



1 

F(A;2) 
1 

1 



+ 



1 



(/c2 -)- p2 



1 



+ 



F(p2) 

1 

F(p2 
F(p2) 



In 



In 



F(A;2) F(p2)y_ 



F(A;2) F(p2 



(143) 

This is our simplest expression for the transverse part. We can then construct the full vertex 
from this using 

4 

T^{p, k-q) = Y. Hp\ k\ q') Lf (p, k;q) + r,(p^ k\ q') Tt{p, k; q) , (144) 

i=l j=2,3,6,8 

from Eqs. (171- fTOl) . This is our final result. Phenomenological studies of strong coupling QED 
with this vertex ansatz are presently underway 5l|, |52| . 



VIII. CONCLUSIONS 



The Schwinger- Dyson equations constitute the field equations of a theory. Being an infi- 
nite set of nested integral equations, they are in general intractable without some form of 
truncation. To date, the only known consistent truncation procedure is perturbation the- 
ory. This satisfies gauge invariance and multiplicative renormalizability order-by-order, and 
the meaning of any truncation is well-defined. In the case of non-perturbative truncations, 
like the rainbow approximation, one has always been unsure as to how much physics has 
been encoded and how much lost. The calculation of dynamical mass generation nicely 
illustrates this. The properties of gauge invariance and multiplicative renormalizability are 
fundamental to our ability to calculate consistently in a gauge theory. It is thus natural 
that any truncation should respect these properties. They ensure not only the elimination 
of overlapping divergences that plague Schwinger-Dyson calculations, but allow all ultravio- 
let divergences to be handled appropriately. Here we have considered the fermion and boson 

50 



propagators in 4-dimensional massless QED. To be able to study these requires an ansatz 
for the full fermion-boson vertex. This interaction involves 11 non-zero components, three 
of which are fixed by the Ward-Green- Takahashi identity in terms of the fermion propagator 
functions. The other eight (transverse) components in principle require knowledge of the 
four, five, six ... point functions. However, very specific projections of this vertex appear in 
the fermion and boson self-energies. We have seen that these projections are strongly con- 
strained by the multiplicative renormalizability of the fermion and boson propagators. At 
its simplest, multiplicative renormalzability is closely related to the ultra-violet behaviour of 
loop integrals. This probes distinct limits for the fermion-boson vertex : one in the fermion 
equation and the other in the boson. In these two limits, the vertex has quite different struc- 
tures. Such behaviour ensures the multiplicative renormalizability of leading logarithms and 
shows that the 2-point Green's functions for both fermion and photon are wholly determined 
by the fermion wavefunction renormalization. This has enabled us to unravel for the first 
time the non-perturbative structure of the full vertex, Eqs. ( 11431 11441 ), at least as far as 
concerns the fermion and photon Schwinger- Dyson equations. 

While the form of the 3-point vertex is determined in three kinematic limits, when k'^,p'^ S> 
g^, when /c^,g^ ^ and when p'^,q'^ ^ k"^, its form at general momenta when all six 
vector structures of massless QED contribute involves free parameters. Imposing the known 
perturbative 0{a) result for the individual vertex components fixes these. This marks a 
significant step in the development of non-perturbative Feynman rules needed for realistic 
calculations in strong QED. There are many steps to go : 



to solve the extended constraints beyond leading logarithmic order and include 
masses |53|, 



to compute the Lamb shift of hydrogen and calculate the properties of positronium to 
asses how well our vertex ansatz automatically sums higher orders in a, 



to explore strong physics with such a complete, unquenched vertex 
existing studies using bare, Ball-Chiu and CP vertices [ll. 



33 



34 



IS 



20 



23 



extending the 



29 



30 



481]. Such calculations are under way and will be reported elsewhere 52] 



31 



Eventually an extension to QCD will be our target. 
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APPENDIX A: PERTURBATIVE r'S 



The vertex coefficients r^'s given below are the mass 
for the massive fermions in general covariant gauge 



ess limit of the exact 0{a) calculation 



4l| 



87rA2' 



^0 



4A2 



1 

+ ln — 



In 



3 



4A2 



^■i^g' + i (e-2) + i 



(Al) 



87rA2' 



8A2 



+ ln- 



+ ln 



{k + pf 



4 

4A2 

(e-2)|, 



(A;2-p2)2 (^_2) 



(A2) 



r^'{p\k\q') 



87rA2 2 



, k' 

+ ln^ 

pZ 



+ ln 



k ■ p{p^ — A;^^ 



4A2 

3 
4A2 



k-pq^ (e-2) 



(p + fc)2 

2(p2 - A;2) 



(e-2) 



+ (e-2) , 



(A3) 



r^\p\k\q') 



87rA2 



{ 



k ■ p Jq + \yi 



]^2p2 



(p2 - A:2) In 



pZ 



(A4) 
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where 



Jo 



2 
A 



k ■ p — A 
p2 



f 



k-p + A 

p2 



In 



In 



k ■ p — A 
k-p + A 



(A5) 



and 



/(x) = Sp{l - x) 



dy- 

1 - y 



(A6) 



APPENDIX B: LIMITS OF Ti'S 



1. Fermion Limit 



In order to take the /c^ ^ ^ p^ hmit of the perturbative transverse vertex coefficients 
namely the Tj functions, Eq. flAmA4l) we need to expand Jo function, Eqs. flA5l IBip . up to 
C(l/F) : 

+ 1 {2{k.pf-^-^{k-pfp' + {k-p)p'- ip^^ 

+ ^ (^-^{k.pf-^-^{k.pfp^ + ^-^{k.pfp'-{k.p)p' + \p'^ 

+ ^ (y(^-;^)^-y(A:-p)V + 10(A:.p)y-^(A:.p)V + (A:.p)p«-^) 

+ 0{l/e) \ . (Bl) 



2. Photon Limit 



In the photon limit, /c^ ^ ^ g^, Jq behaves like : 



^0 



(p2 - A;2) 



A; 



2\2 



p'i 



13 - ky 
18 / 



(B2) 
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APPENDIX C: EFFECTIVE r'S 



The connection between the effective and real Tj functions are given below and the detail of 
this procedure can be found elsewhere [50| : 



2 _ 



where 



1 ''\^'_!!^^riUUp\k\q') \-A^-lqH.p] 



f{k\v') Jo 
1 



{r,%ff{p',k^) = y^T^ r d^'-^{r,^)neaiip',k^q') {k ■ p} , 
Jeyi^ jP ) Jo H 



.2 ^2n _ 



f{k\p 
k{k\p' 



8 ^ 
IT k"^ 

4 F ■ 



(CI) 



(C2) 
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